Controlling entanglement by direct quantum feedback 
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We discuss the generation of entanglement between electronic states of two atoms in a cavity 
using direct quantum feedback schemes. We compare the effects of different control Hamiltonians 
and detection processes in the performance of entanglement production and show that the quantum- 
jump-based feedback proposed by us in Phys. Rev. A 76 010301(R) (2007) can protect highly 
entangled states against decoherence. We provide analytical results that explain the robustness of 
jump feedback, and also analyse the perspectives of experimental implementation by scrutinising 
the effects of imperfections and approximations in our model. 
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I. INTRODUCTION 

Recent experimental advances have enabled individ- 
ual systems to be monitored and manipulated at the 
quantum level in real time [llj S S 0, S E|- As new 
advances continue to emerge, and the development of 
control strategies for quantum systems becomes essen- 
tial, quantum technology inevitably approaches the well- 
developed classical control theory, borrowing its concepts 
and extending them to the quantum realm. An impor- 
tant example is feedback control, which consists of the 
manipulation of the system according to the informa- 
tion acquired through measurement. Feedback control 
has been applied to quantum systems 0, [i, [3, [lO|, in- 
cluding its irnplementation in a variety of experimental 
setups [l,[il|,[il,l3- 

Relying on the ability to produce specific states and 
perform controlled operations on them, quantum infor- 
mation is undoubtedly an area that would benefit from 
further advances in quantum control. Of particular in- 
terest in this context is the preparation of entangled 
states, indispensable for quantum information process- 
ing. Notwithstanding the considerable number of ex- 
periments carried out on the generation of entangle- 
ment [ll,[ll,[ll,[I3,IIl|, and the effort to protect the sys- 
tem against undesirable imperfections and interactions, 
entanglement decay due to uncontrolled coupling with 
the environment remains a major problem yet to over- 
come [l^ [20, |21|. In this scenario, quantum feedback 
emerges as a possible route to develop strategies to cir- 
cumvent entanglement deterioration. 

In fact, quantum feedback control has been recently 
used to improve the creation of steady state entangle- 
ment in open quantum systems. A Markovian, or di- 
rect, feedback [^, Q was used in ]2^, [23| to show that 
states with entanglement corresponding to one third of 
the maximum possible value can be produced. Maxi- 
mally entangled states could be achieved [2J] in an ide- 
alised situation, by the use of Bayesian [l^l, or state esti- 
mation, feedback. This improvement, however, comes at 
the cost of an increasing experimental complexity due to 



the need of a real time estimation of the quantum state in 
the later method, as compared to the simple feedback di- 
rectly proportional to the measurement signal proposed 
by the former strategy. 

Despite being simpler to implement, direct feedback 
still exhibits a multitude of possibilities due to the ar- 
bitrariness of choices for the control Hamiltonian and 
measurement schemes. In a recent Rapid Communica- 
tion [25|, we have shown that an appropriate selection 
of the feedback Hamilton and detection strategy leads to 
the robust production of highly entangled states of two 
atoms in a cavity. In this contribution, we will further 
explore the richness of feedback strategies by comparing 
the jump-based scheme proposed in 12511 with a strategy 
based on homodyne measurements |23|, taking into ac- 
count the unavoidable imperfections that may occur in 
realistic experimental implementations. 

The paper is organized as follows. In the first part of 
Section |TT] we introduce the model of the system and its 
equation in the absence of feedback. In the second part 
we describe the introduction of feedback and the corre- 
sponding changes in the master equation for two different 
choices of measurement strategy, namely, photodetection 
and homodyning. Entanglement dynamics is examined 
in Section lllli and the role of feedback on entanglement 
generation is analysed. Robustness issues are discussed 
in Section [TVl where the effects of detection inefficiencies 
and spontaneous emission on steady-state entanglement 
are considered. Section FVl discusses different experimen- 
tal regimes and explores the non-adiabatic limit of the 
model. Section IVII concludes with a discussion of our 
main results. 



II. THE MODEL 

The system consists of a pair of two-level atoms cou- 
pled resonantly to a single cavity mode, with equal cou- 
pling strength g, and simultaneously driven by a laser 
field with Rabi frequency O. The cavity mode is damped 
with decay rate k and the atoms can spontaneously decay 
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FIG. 1: Schematic view of the model. The system consists of 
a pair of two-level atoms coupled to cavity, which is driven by 
a laser field and damped. Conditioned on the measurement 
of the output of the leaky cavity, a Hamiltonian is applied to 
the atoms, completing the feedback scheme. 



with rates 71 and 72. Feedback is applied conditioned 
on the measurement of photons leaving the cavity, as 
schematically shown in Fig. [1] In the absence of feed- 
back, the master equation describing the system is given 
by 

p ^ -m [{J+ + J^),p] - ig [{J+a + J-a''),p] 



-KV[a]p + Y,rM<^^]P■ (1) 



The first and second terms represent the Hamiltonian 
evolution induced by the laser driving and atom-cavity 
coupling, respectively. The superoperator 
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I?[c]p = cpc* — - (^c^cp + pc'c) 



(2) 



describes cavity and atomic decays, given in terms of a, 
the annihilation operator of photons in the cavity, and 
ai = \gi){ei\, the lowering operator for the i-th two level 
atom, respectively. The angular momentum operators 
are defined as 



J- 



0-1 +0-2, 
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(3) 

(4) 



with Gi 



\ei){gi\ the raising operator for the atomic 



electronic levels. 

In the limit where the cavity decay rate k is much 
larger than the other relevant frequencies of the problem, 
the cavity mode can be adiabatically eliminated and one 
obtains a master equation just for the atomic degrees of 
freedom [23| 

p = ~in[{j+ + j_),p] + rp[j_]p + ^7,i?Hp, (5) 

i 

where T — g^ j k. is the effective collective decay rate. Fur- 
thermore, the Dicke model [2y| can be recovered under 
the assumption that the collective decay rate is much 
larger than the spontaneous emission rates, F ^ 71, 72, 



Later we shall relax the approximations used to ob- 
tain Eqs. (O and (O: The non-adiabatic regime is im- 
portant when considering experimental situations where 
high quality cavities are used, so it will be discussed in 
Section |Vl Section IIVI will show that even small sponta- 
neous emission rates can have a drastic deleterious effect 
on the final amount of entanglement produced in this sys- 
tem. For the moment, Eq. ^^ will be the starting point 
to introduce the description of the measurement scheme 
and the feedback mechanism. 

Here we will consider the direct, or Markovian, feed- 
back introduced by Wiseman and Milburn [^, Q , where 
the control Hamiltonian is proportional to the measure- 
ment signal. This control mechanism has the advan- 
tage of being simple to apply in practice, since it avoids 
the challenge of real time state estimation required in 
Bayesian, or state-based, feedback [T^. The idea is de- 
picted in Fig. [T] the cavity output is measured by a de- 
tection scheme M whose signal I{t) provides the input to 
the application of the control Hamiltonian Tffb = /(t) F . 
Here, we will consider the measurement stage M to be ei- 
ther a homodyne or a direct photodetection of the output 
field. 

In the homodyne-based scheme, the detector registers 
a continuous photocurrent, and the feedback Hamilto- 
nian is constantly applied to the system. Conversely, in 
the photocounting-based strategy, the absence of signal 
predominates and the control is only triggered after a 
detection click, i.e. a quantum jump, occurs. This is re- 
flected in the different forms of the equations representing 
the dynamics of a feedback-controlled system under ei- 
ther of the measurement schemes. In the homodyne case, 
Eq. © becomes [1, ij 



p^ ChP^~-^[iJ^ 



j-),p] + rv[j^]p + ^v[F]p 

-^[F,-iJ^p + ipJ+].i7) 



Cp = -iVL [{J+ + J_), p] + FP[J_]p. 



(6) 



This equation was used by Wang, Wiseman and Milburn 
(WWM) [111 to show that the amount of steady state 
entanglement in this atom-cavity model can be enhanced 
above that generated by the uncontrolled dynamics [23] . 
More recently, the same equation was used to explore the 
effect of different feedback Hamiltonians on entanglement 
generation [28J. 

In a measurement scenario based on photodetec- 
tions, the unconditioned dynamics, including feedback, 
reads 

p = CpP = -'-n [(J+ + J_),p] + TVpJ^p. (8) 

The manifestation of the abrupt character of the jump 
feedback is clear when one writes the last term of Eq. ([5]) 
explicitly: V[UJJ\ = UJ-pJ+W - {J+J-p + pJ+J-)/2. 
The unitary transformation U — exjp[—iF6t/h], repre- 
senting the finite amount of evolution imposed by the 
control Hamiltonian on the system, only acts immedi- 
ately after a detection event, which is described by the 



first term of the superoperator V. Note that this is equiv- 
alent to replacing the jump operator J_ in the master 
equation by UJ-. In this way, one can engineer the 
reservoir dynamics [2^, [30] via the transformation U to 
generate highly entangled states. 

From the differences between Eqs. ^ and ^ it is evi- 
dent that the resulting dynamics should strongly depend 
on the choice of measurement strategy. Here, we will ex- 
plore this fact, together with the possibility of changing 
the feedback Hamiltonian to show how they affect the dy- 
namics of the system and, in particular, its asymptotic 
entanglement. 



III. STEADY STATE ENTANGLEMENT AND 
FEEDBACK 




A. Dynamics without feedback 

The steady state solutions for the Dicke model, Eq. ([6]), 
were obtained already many years ago [3lLl32l|. However, 
only recently their entanglement properties were brought 
to attention [27| and explored in the context of quantum 
information. 

The first important feature of Eq. ^ is the fact that 
it is symmetric with respect to exchange of the atoms. 
This suggests that, instead of using the two qubit basis 
{\99)j l^e), \eg), |ee)}, one should use angular momen- 
tum states, and analyse the system in terms of the sym- 
metric (j = 1) 
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1.9.9), |2) 



l.9e) + \eg) 



V2 



, |3) 



and anti-symmetric (j ~ 0) 



|4) 



|ge) - |e.9) 
V2 



(9) 



(10) 



subspaces. 

The anti-symmetric subspace is a decoherence-free 
subspace [3J, |33| and the state |4) is therefore a steady 
state solution of Eq. ©. This is an interesting case, de- 
spite its trivial dynamics, since the asymptotic state in 
this subspace is a pure, maximally entangled, one. In 
fact, this situation was explored in a recent proposal for 
producing Werner states in a system of atoms inside a 
cavity [35| and in a probabilistic scheme to generate the 
singlet state via quantum-jump detection [36|. In terms 
of dynamics, however, the symmetric case is more inter- 
esting: while in the anti-symmetric subspace an initially 
prepared Bell state |4) does not evolve at all, entangle- 
ment can be dynamically generated from any symmetri- 
cal initial condition, even from initially separable states. 
However, even for optimal parameters, the amount of 
entanglement in this case is only about 10% of the Bell 
state's value |27| . 

For a general asymmetric initial condition the situa- 
tion gets more complicated as both symmetric and anti- 
symmetric components are present. This is exemplified 



FIG. 2: Steady state concurrence as a function of the ratio 
Sl/r for initial the states \ge) (solid line) and \gg) (dashed 
line). In the anti-symmetric subspace there is no evolution 
and the system is always in the maximally entangled state 
4) (not shown), while for the symmetric part the maximum 
entanglement is reached for Q/T ~ 0.38 [23]. The steady state 
for the asymmetric initial condition strongly depends on the 
interplay between symmetric and anti-symmetric solutions. 



in Fig. [21 where the steady state entanglement, measured 
by the concurrence [37], as a function of the ratio fl/T is 
shown for an asymmetric {\ge), solid line) and a symmet- 
ric {\gg), dashed line) initial states. For zero driving, the 
symmetric part of the asymmetric state evolves towards 
the ground state {\gg)) and the entanglement is totally 
due to the anti-symmetric component (]4)). Tuning the 
parameters to obtain the maximum entanglement for the 
symmetric component (fi/F w 0.38) leads to a smaller 
amount of entanglement in the asymmetric case. This 
happens because in the latter case the solutions arising 
from different subspaces interfere, eventually resulting in 
a less entangled state. 



B. Dynamics with feedback 

Inclusion of Markovian feedback can significantly in- 
crease the amount of asymptotic entanglement as com- 
pared to the situation described in Fig. [51 This was 
first shown in [23] for a homodyne-based feedback where 
the steady state entanglement was approximately 3 times 
larger than in the non-controlled case. Also in the case of 
homodyne detection, a steady-state entanglement of 80 
percent of the maximal possible value was obtained 12811 
using the local asymmetric feedback law proposed in ^^ . 
In this section we shall show how to improve these results 
by designing different feedback schemes. 







FIG. 3: (Color online) Steady state concurrence as a function 
of driving and feedback strengths (left) for homodyne (a) and 
photo-detection feedback (c). The homodyne case reproduces 
the result of I23| with a maximum concurrence c ~ 0.31. Us- 
ing a jump feedback strategy, the maximum entanglement 
increases to c « 0.49. On the right, the absolute value of the 
density matrix elements (in the angular momentum basis) for 
the parameters corresponding to the maximum entanglement 
in the left panel. 



1. The role of measurement: homodyne vs. 
photodetection-based feedback 

We begin our analysis recalling the WWM results [2j| . 
Their scheme fits the general scenario depicted in Fig. [1] 
when the measurement stage M is a homodyne detection 
and the feedback Hamiltonian is F = XJ^ = X{J- + J+). 
In fact, the form of the control Hamiltonian coincides 
with the one appearing in Eq. ^ and in 23] was pro- 
posed to be implemented via a modulation of the field 
driving the cavity. With this choice of control, the final 
equation, Eq. ([7]), retains the symmetry properties with 
respect to exchange of atoms. Assuming a symmetric 
initial condition, the system will remain in the subspace 
given by Eq. ^ and an analytical solution for the steady 
state can be found [23| ■ Figure [3^ shows the entangle- 
ment of these solutions as a function of the driving and 
feedback strengths, with a maximum concurrence c ~ 0.3 
obtained for n/T « ±0.4 and A/F w -0.8. 

Now, we change the measurement process from homo- 
dyne to photo-detection and investigate the effects of this 
replacement on entanglement generation. For the mo- 
ment, we keep the same form for the control Hamilto- 
nian in order to focus on the changes induced solely by 
the alteration of the detection scheme. Note that, as 
mentioned before, the feedback now is applied only when 
a detection occurs and is described by the action of the 
operator U = exp[—iFSt/h] = exp[— lAJ^:]. The asymp- 
totic solution in the symmetric subspace can be calcu- 



lated from Eq. ([8]) and the corresponding entanglement 
is shown in Fig. [3J;. Now, the maximum concurrence, 
c « 0.49, occurs at fl/T w -0.0023 and A/F w ± - 1.57, 
exceeding the value obtained via homodyne-based feed- 
back. 

The sharp differences between the two feedback 
schemes are evidenced by the contrasting plots of Fig. [3] 
(a) and (c) and by the form of the steady states (at the 
peaks) depicted in Fig. [3](b) and (d), where the absolute 
values of the density matrix elements in the angular mo- 
mentum basis, Eqs. ^ and PH)) . are shown. In the jump 
feedback case, the final state at the peak concurrence is 
close to a mixture between |3) and the Bell state |2), 
the latter being the state responsible for the entangle- 
ment. In the homodyne case, the abundance of different 
off-diagonal elements suggests a more complicated struc- 
ture. Understanding these differences will be important 
for the analysis of spontaneous emission effects in Sec- 
tion EVl 



2. The role of the control Hamiltonian 

A change in the feedback Hamiltonian to improve the 
control introduces an enormous range of new possibili- 
ties, even when considering the limitations imposed by 
constraints on experimental feasibility. In this section, 
however, we will restrict the discussion to the case of 
a local feedback, i.e. when a control Hamiltonian act- 
ing on just one of the atoms is employed, as proposed 
in [25|. The reason for this choice is that a local feed- 
back breaks the symmetry of the system, making it possi- 
ble to explore the interplay between symmetric and anti- 
symmetric subspaces from a dynamical point of view. 
Under the symmetric control described previously, all 
symmetric initial states lead to the same final station- 
ary solution, while the anti-symmetric state is station- 
ary itself. There are, however, infinitely many stationary 
states mixing both subspaces, which are unambiguously 
determined by the choice of initial conditions. The first 
important point to address is then how the relaxation 
properties of the feedback master equations change un- 
der the new asymmetric control law. 

Our discussion here will be based on the coherence- 
vector formalism as described by Lendi in [3^. In this 
approach, the original master equation is transformed 
into a linear system of equations in real space, i.e. 
p{t) = Cp{t) -^ v{t) = Gv{t) + k, and the stationarity 
properties will be encoded in the matrix G. A well known 
example of such a transformation is the description of 
two level systems in terms of the Bloch equations. For 
our purposes it suffices to note that if det(G) ^ then 
the system admits a unique stationary solution given by 
Vss = —G^^k. If dct(G') = then two situations may 
arise: either the equation Gv{t) + k = has no solution 
at all, or it has infinitely many, determined by the initial 
condition. 

The matrix G can be straightforwardly calculated from 



the feedback master equations and it is obvious, from our 
previous discussion on the symmetric control J^, that, in 
that case, G is singular since there were more than one 
stationary solution. However, for a general local control 
of the form U = Ui®\, the matrix is non-singular (except 
for the trivial cases oi Ui — \ and f7 = 0) and there is 
a unique solution. Moreover, note that, in the case of 
jump-based feedback, the anti-symmetric Bell state |4) 
remains stationary for any choice of U and therefore is 
the only steady state of the system. 

In the case of homodyne feedback, the steady state 
solution is also unique, but depends on the form of Ui 
and on the parameters A and Q.. Figure [4] illustrates the 
effect of local control on the steady state entanglement 
for the particular choice 



F^Xa,. 



(11) 



with ax = (Jx+ai , for homodyne (a) and photo-detection 
(c) based feedback. In the homodyne case, there is a 
substantial increase in the asymptotic entanglement as 
compared to the Jx control of Fig. [31 The final state with 
the largest amount of entanglement corresponds now to 
a combination of the anti-symmetric component |4) with 
the symmetric ones, |1) and |2) (see Fig. Wp), with c = 
0.81 for A/F = ±0.01 and n/T = ±0.07. As discussed in 
the previous paragraph, for a photo-detection feedback 
based on a local control law, a maximally entangled Bell 
state is generated from any initial condition for all non- 
trivial parameters, as shown by the plateaux in Fig. |4l; 
and the density matrix elements in Fig. |4li. 

Since this steady state can also be obtained directly 
from the non-controlled system, one could, at first sight, 
question the relevance of the use of feedback in this case. 
Indeed, if one considers the model described by Eq. ([6|) 
without any kind of imperfection, then the jump feedback 
with local control would represent no advantage over the 
non-controlled case. Nonetheless, imperfections are in- 
herent to any real physical system, and should be con- 
sidered. Without control, for example, the asymmetric 
Bell state has to be produced beforehand to be unaf- 
fected by the dynamics, and any symmetric component 
introduced by a non-ideal preparation would spoil the fi- 
nal entanglement (see Section Fill A[) . Conversely, with 
feedback the system naturally evolves to the pure entan- 
gled state |4) for any initial condition. However, besides 
non-ideal initial state preparation, other sorts of imper- 
fection as, for example, detection inefhciency or errors in 
the production of the feedback Hamiltonian, may arise 
from the feedback implementation, and will investigated 
in the next section. 



IV. ROBUSTNESS OF CONTROL 

A. Spontaneous emission effects 

All our previous discussions were based on the appli- 
cation of different kinds of feedback control on the model 




FIG. 4: (Color online) Steady state concurrence for a homo- 
dyne (a) and photo-detection (c) feedback with local control 
(Ta; . There is a huge improvement as compared to the situation 
depicted in Fig. [31 For a local control, a pure anti-symmetric 
Bell state is obtained for all parameters (except in the cases 
with zero driving SI = or no feedback A = 2-nm), and all 
possible forms of U\. On the right, the absolute value of the 
density matrix elements for the peaks on the left plots. 



described by Eq. ^ , where spontaneous emission effects 
were neglected. However, even if all other sources of im- 
perfection were surmounted, spontaneous emission would 
still be the fundamental limiting factor for the existence 
of entanglement in a system of atomic qubits. Conse- 
quently, the ultimate goal for the feedback schemes in- 
vestigated here would be the production of steady state 
entanglement under a model that includes this effect. 

The starting point is therefore Eq. ([5]). The new feed- 
back equations are similar to Eqs. ([7]) and ([5]), differing 
only by the addition of the last two terms of Eq. ([5]). 
The first thing to note is that those terms break the 
decoherence-free condition for the anti-symmetric sub- 
space and the advantages of using feedback will become 
even more evident. Analytical solutions for the steady 
states can still be found but, being generally cumber- 
some, will not be presented here in their complete form. 
Instead, we will analyse the regime where spontaneous 
emission effects can be considered as a perturbation to 
the system, and focus on the qualitative understanding 
of how these effects influence the feedback schemes dis- 
cussed in the previous section. 

Let us consider first the scheme with symmetric con- 
trol where, in the absence of atomic decay, symmetric and 
anti-symmetric subspaces were decoupled and the system 
presented infinitely many stationary solutions. Sponta- 
neous emission, however, introduces a coupling between 
the subspaces as it pushes the system towards its ground 
state. During this process, the anti-symmetric compo- 
nent 1 4) is produced, breaking the symmetry of the sys- 



tern. This has a huge impact on the new steady state of 
the system, which is now unique and, for most of the pa- 
rameters, has no resemblance to the corresponding state 
for 7 = 0. The corresponding steady state entanglement, 
shown in Fig. [5};, is drastically reduced as compared to 
the situation of Fig. [8)3, despite the small value of the 
atomic decay rate. 

The scenario changes drastically when the local feed- 
back Hamiltonian pT|) is considered. In this case, even 
without spontaneous decay, the system could explore the 
whole Hilbert space as the original symmetry was al- 
ready broken by the control Hamiltonian. Moreover, the 
steady states, for both homodyne and jump-based feed- 
back, were unique and highly entangled. The results for 
the steady state entanglement for 7 = 0.01 F are shown in 
Figs. [5)3 and[5jl for homodyne and jump feedback, respec- 
tively. Differently form the Jx control case, spontaneous 
emission has indeed only a perturbative effect, reducing 
slightly the maximum value of steady state entanglement 
(compare with Fig. 2]). The figure also indicates that, for 
the local jump control, the steady state remains close to 
the anti-symmetric Bell state for most values of A and 
fi, preserving the plateau structure shown in Fig.UJ) (see 
Fig. [SJi). In fact, this can be shown directly from the 
analytical solution including atomic decay. The full solu- 
tion is complicated but, expanding in powers of 7/r and 
keeping only the lowest order term, one obtains 




Jump 



FIG. 5: (Color online) Spontaneous emission effects on the 
steady state concurrence for the cases corresponding to Fig. |3] 
(left, Jx conrol) and Fig. U (right, g^; control). Even for a 
small decay rate (7/r = 0.01), the effects for a J^ control 
are pronounced for both homodyne (top) and jump (bottom) 
detection schemes. In the local control case (right) the shape 
seen in Fig. |4]is left basically unaltered with only a small de- 
crease in the maximum value of concurrence. Remarkably, for 
the local control Eq. (|11[) with photo-detection (d), a highly 
entangled state is stabilised for almost all parameter space. 



B. Detection inefficiencies effects 



P44 « l-0(7/F), 

Py « 0(7/F), z,j^4. 



(12) 



The final state is therefore close to the anti-symmetric 
Bell state. The first order terms depend on the parame- 
ters ri and A, and this dependence is responsible for the 
small deviations from the plateau structure of Fig. 0):. 

The major reason for this robustness of a local control 
strategy using both detection schemes lies in the unique- 
ness of the steady state without decay. For any initial 
condition, the system is driven to the target state and, 
therefore, as soon as spontaneous decay forces the system 
away from this state, feedback dynamics counteract this 
tendency. This competition between feedback and decay 
dynamics determines how far the perturbed stationary 
state will be from the original one. Consequently, as soon 
as the ratio between atomic and collective decay rates re- 
main small, the final state will remain highly entangled. 

These results are not restricted to a feedback Hamilto- 
nian of the form Eq. (|lip . In fact, similar steady state en- 
tanglement arc obtained other local control laws, though 
less entanglement can also be observed for particular re- 
gions of parameters. Therefore, there is still room for 
improvements in the results shown in Fig. [5] as one can, 
in principle, optimise over the forms of U to get even 
larger final entanglement. 



Since feedback relies on the manipulation of the system 
based on information gained by a measurement, evalua- 
tion of the effects of inefficiencies in the detection pro- 
cess is important. Finite efficiency can be introduced 
both for photo-detection [39| and homodyne measure- 
ments |40|, and this has been explicitly considered for 
homodyne feedback schemes [3, |4l|, |42|. From now on, 
we will consider only the jump-based feedback, since it is 
the scheme with best performance for entanglement gen- 
eration. In this case, the extension of Eq. ([8]) to allow 
for a non-unit-cfficicncy detection can be done by identi- 
fying two distinct situations when a jump occurs: in the 
first the detector clicks and the feedback transformation 
U acts in the system, in the second the detector fails to 
click and no control is applied. The corresponding equa- 
tion reads 



--n [{J+ + J_), p] + Tr^VpJ^p 

h 

+F(l-7?)I?[J_]p + ^7,PHp. 



(13) 



When the detector efficiency rj is zero, no information 
is extracted from the nreasurement and the equation re- 
duces to the equation without feedback, Eq. ([5]). Ev- 
idently, for an unity efficiency detector Eq. ([5]) is re- 
gained, and, for a local control, a maximally entangled 
steady state is reached. In the intermediate case where 
< 77 < 1, one would expect that imperfect knowledge 
gain should lead to inefficient control. Note however that. 




FIG. 6: (Color online) Maximum steady state concurrence as 
a function of the ratio 7/r and detection efficiency rj. For 
small 7/r the system is almost insensitive to detection in- 
efficiencies, decaying abruptly when r/ approaches zero (for 
7/r = entanglement is maximum unless rj = 0). 



neglecting spontaneous emission, the anti-symmetric Bell 
state is a steady state of both Eqs. ^ and ^, and 
one can show, proceeding exactly as in the unit-efRciency 
case, that this also holds true for Eq. p3)) for any 77 > 0. 
This is a peculiar situation of this dynamics: while the 
feedback term tries to move the system to the state |4), 
the one corresponding to a missed click does not affect 
this anti-symmetric component, impinging only a delay 
in the time taken to reach the steady state [2^| . 

This is not true if atomic decay is taken into account 
as undetected events imply missed opportunities to ap- 
ply the feedback, hence an effectively weaker feedback 
as compared to spontaneous emission. Mathematically 
this is clear from Eq. P^ as the feedback term is scaled 
by rj. The effect of detector inefficiencies and sponta- 
neous emission together are summarised in Fig. [6] where 
the maximum steady state concurrence is plotted as a 
function of rj and 7/r. The most noticeable feature is 
the weak dependence of the entanglement on the detec- 
tion efficiency in the limit 7 ^ F (unless 77 is close to 
zero) that confirms the robustness of a local jump-based 
feedback in this regime [25|. For a ratio 7/r = 0.002, for 
example, one still has concurrence above 0.9 for detection 
efficiency as low as ry = 0.1 



V. ANALYSIS OF THE ADIABATIC 
APPROXIMATION 

The analysis from the previous sections shows that en- 
tanglement production can be vastly improved by the 
use of feedback strategies and that, in the best case sce- 
nario, highly entangled states can be protected against 
decoherence using a quantum-jump-based feedback with 



a local control Hamiltonian. The robustness of the feed- 
back scheme depends on a single quantity, the ratio 
C — r/7_= 5^/k7, also known as the cooperativity pa- 
rameter [43| . The higher the cooperativity, the better the 
performance of the feedback procedure: not only does the 
steady state entanglement obtained increase, but also the 
sensitivity of the scheme to detection inefficiencies is re- 
duced. 

Interestingly, large cooperativity, which corresponds to 
the strong coupling limit, is the regime used by many of 
the current cavity QED experiments in the optical do- 
main p, |4J, |45|. This is achieved either by increasing 
the atom-cavity coupling or decreasing the cavity decay 
rate, since the spontaneous emission rate is fixed by the 
choice of atomic species. The use of high quality cavities, 
however, affects the feedback model based on Eq. (O, 
which is obtained in the bad cavity limit after adiabatic 
eliminating the cavity mode. To have a better under- 
standing of how our model would work within current 
experimental conditions and also to determine the most 
favourable region of parameters for future experiments, 
in this section we will relax the adiabatic condition and 
analyse the performance of local jump feedback strategy 
in this situation. 



A. The non-adiabatic model 

Equation ^ describes the whole system, including 
cavity mode dynamics, without control. Introduction of 
photo-detection feedback follows exactly the same rea- 
soning used to obtain Eq. ([8]). Whenever a photon is de- 
tected, the control Hamiltonian is applied and the system 
undergoes an evolution given by the operator U. This 
operator now enters in the term that describes the mon- 
itoring of the cavity output, i.e. I?[a]p. Addition of de- 
tection inefficiency also follows directly from our previous 
discussion and the full master equation for the system is 

p = -ifl [{J+ + J_),p] - ig [{J+a + J-a^),p\ 

+riKD[Ua]p + (1 - ri)KV[a]p + ^ l^V[a^]p, (14) 

i 

which transforms back to Eq. (|13p if an adiabatic elimi- 
nation is performed. 

To start with, instead of solving the full master 
equation. E g. (fM]) . we will use a quantum trajectory 
method [4^ |43| to follow the dynamics conditioned on 
the measurement results. This approach, which provides 
a useful picture of single runs of an actual experiment, 
combines continuous evolution under an effective Hamil- 
tonian and the random occurrence of quantum jumps. In 
our problem there are four possible jumps, corresponding 
to the four decoherence terms in Eq. p^ : a detected pho- 
ton leaving the cavity that triggers the feedback (third 
term), an undetected photon (fourth term), and sponta- 
neous emission from the atoms (two jumps represented 
by the sum in the fifth term). For simulation purposes, 



we assume that the environment is perfectly monitored 
so that the state remains pure throughout the evolution. 
Obviously this is not true for undetected photons, but 
we can take a different perspective and consider feed- 
back inefficiencies as detected photons where the feed- 
back Hamiltonian fails to apply. Note that for the aver- 
age behavior given by Eq. (I14p those points of view are 
equivalent. 

Figure [7] shows the concurrence (top panels) and the 
average number of photons in the cavity (bottom panels) 
for single stochastic trajectories as a function of time. 
The cooperativity parameter is set to C = 100 (as in 
Fig. [5]) , the detection efhciency is ry = 1 , and we assume 
that the atoms are initially in the ground state and the 
cavity in the vacuum. Figure [7^ corresponds to a situ- 
ation close to the adiabatic regime with k = 4OO7 and 
g = 2OO7, while g = 4O7 and k = I67 in Fig. Uh- AH 
other parameters are selected to give the maximum con- 
currence for the chosen g and k. The average behavior, 
given by the solution of Eq. [Ml is shown by the dashed 
lines for comparison. 

For the case of Fig. [7^, dynamics with feedback drives 
the system to the target state in a short time. In this 
particular realization, two detected photons, and con- 
sequently two control pulses, at i « 0.15 are enough 
to correct the trajectory to state |4) with no photons 
in the cavity. The absence of photons from the cavity 
output for certain periods of time is therefore an in- 
dication that the atoms are in the anti-symmetric Bell 
state. This state can only be perturbed by the oc- 
currence of spontaneous emission jump (in this case at 
-ft « 0.62,1.15,2.3,2.8,2.9, and3.5), which destroys en- 
tanglement. As soon as the system gets away from the 
dark state, cavity is again populated and control pulses 
applied when photons are detected, reestablishing the 
maximally entangled state. In the adiabatic regime, since 
K 3> 7, this happens quickly during short windows of 
control jumps (we have 7 of them in Fig. [7^). As a re- 
sult, the system spends most of the time in the entangled 
state and, on average, entanglement will be high. This 
is another way of understanding the perturbative nature 
of the spontaneous emission in the adiabatic limit: the 
larger the cooperativity, the larger the amount of time 
spent in the target state, hence higher entanglement. 

However, this picture is only valid in the adiabatic limit 
where cooperativity is the only important parameter. 
Fig. [7|3 shows the dynamics outside the adiabatic regime 
but for the same cooperativity. Now, with a smaller k, 
the windows of control are larger and the system takes 
longer to reach the target state. Moreover, the ratio be- 
tween the probabilities of having spontaneous emission 
and control events increases. Consequently, the system 
spends less time in the state |4), resulting in a less entan- 
gled steady state (c ~ 0.6). This may be explained by 
the reliance of the control scheme on information trans- 
mitted by the photons that leak from the cavity. In the 
same way that low detection efficiency limits the eventual 
peak concurrence because an opportunity for control is 
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FIG. 7: Concurrence (top panels) and mean cavity photon 
number (bottom panels) as a function of time for single quan- 
tum trajectories. In (a) the system is in the adiabatic regime 
with K = 4OO7, g = 2OO7, n = 4O7 and A = -I87. In (b), 
effects of non-adiabatic dynamics get more evident (k — I67, 
g = 4O7, n = 20.87 and A = -IO.O5287). In both plots, 
the cooperativity is C = 100, and the concurrence for the 
time-evolved density matrix is shown for comparison (dashed 
lines) 



lost, if the cavity decay rate is very low, the photons that 
are emitted into the cavity have a higher probability of 
being reabsorbed by the atoms rather than escaping and 
triggering a control pulse. 

One can therefore conclude that the major aspect re- 
sponsible for the robustness of the local jump-based feed- 
back scheme is the balance between control and decay 
dynamics. In the adiabatic case, the feedback time scale 
is dictated by the effective rate F and the robustness 
depends only on the cooperativity r/7. In the non- 
adiabatic case, feedback rate is given by the cavity decay 
K and the feedback performance will crucially depend on 
the ratio K/7. This indicates that the adiabatic regime 
is the best one to achieve a robust generation of highly 
entangled states. 

Note, however, that the ratio K/7 is not the only im- 
portant parameter in the non-adiabatic case. The laser 




FIG. 8: (Color online) Steady state concurrence for jump- 
based feedback as a function of driving and feedback strengths 
for g — 407,K = I67, and rj = 0.5. The maximum concur- 
rence (c ~ 0.5) occurs for specific values of the parameters 
in contrast with the flat structure observed in the adiabatic 
regime (see Fig. [SJl). 



can excite the atoms, which can then emit photons in the 
cavity mode. These photons can finally escape the cav- 
ity, triggering the feedback. The dynamics depends, in 
an intricate way, on how the rates at which those process 
occur are related, contributing significantly for the be- 
havior of the steady state entanglement. Figure [5] shows 
the stationary concurrence as a function of 51 and A for 
the same values of g and k as in Fig. [71d, but for a smaller 
detection efficiency 77 = 0.5. The simple plateau feature 
from the adiabatic case (Fig.^Ji) disappears, giving place 
to a structure that highlights the strong dependence of 
the final entanglement on the different parameters of the 
problem. 



schemes to the success of the control, showing that a 
photo-detection feedback strategy outperforms the one 
based on homodyne measurements. The effects of dif- 
ferent feedback Hamiltonians were also analysed: local 
control schemes, i.e feedback applied to only one of the 
atoms, produced larger amounts of steady state entan- 
glement than collective interaction. We identified the 
strategy of jump-based feedback with local control [23 
to be the best for robust preparation of highly entangled 
states. For this scheme, we extended the numerical anal- 
ysis given in [2^, and presented a justification for the ro- 
bustness against spontaneous emission and detection in- 
efficiencies based on analytical considerations. Motivated 
by the parameters used in most of the recent optical cav- 
ity QED experiments, we also investigated the role of the 
adiabatic elimination on our feedback scheme and showed 
that the adiabatic regime gives the best performance in 
terms of entanglement generation and robustness. 

An important issue is whether our feedback could be 
realized experimentally. The setup of two atoms equally 
coupled to a cavity mode with possibility of individ- 
ual addressing has already been demonstrated in }48j . 
The remaining constraints for an efficient feedback pro- 
cedure are: i) large cooperativity parameter (F ^ 7) 
and ii) high detection rate as compared to the decay rate 
(k ^ 7). Although the later condition can be certainly 
achieved in experiments with low quality cavities, and the 
former has been obtained in a variety of recent experi- 
ments [4^, [so, [Sll, I54I , the difficulty relies on combining 
all those ingredients in a single experimental setup. The 
usual strategy to fulfill (i) is to improve the quality of 
the cavity, which would be detrimental to condition (ii). 
The solution would be to design a cavity with high trans- 
missivity (ii), yet with a large coupling strength g, which 
would involve a compromise between small mode volume 
and the requirement of individual addressing of atoms. 



VI. CONCLUSIONS 
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